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Abstract 

We examine the statistical mechanics of a 1-dimensional gas of both adjoint and fun- 
damental representation quarks which interact with each other through 1+ 1-dimensional 
U{N) gauge fields. Using large-A^ expansion we show that, when the density of funda- 
mental quarks is small, there is a first order phase transition at a critical temperature 
and adjoint quark density which can be interpreted as deconfinement. When the funda- 
mental quark density is comparable to the adjoint quark density, the phase transition 
becomes a third order one. We formulate a way to distinguish the phases by considering 
the expectation values of high winding number Polyakov loop operators. 



1 Introduction 



The classical Coulomb gas is an important model in statistical mechanics. It is exactly 
solvable in one dimension. In two dimensions it exhibits the Berezinsky-Kosterlitz-Thouless 
phase transition which is the prototype of all phase transitions in two-dimensional systems 
which have U{1) symmetry. In this paper, we shall discuss a generalization of the classical 
Coulomb gas to a system of quarks which interact with each other through non-Abelian 
electric fields. This model is known to be exactly solvable in some special cases, for SU{2) 
gauge group and fundamental representation quarks in one dimension and for SU{N) 
gauge group in the large- limit with adjoint representation quarks in one dimension |0, ^. 
It can also be formulated on the lattice and solved with adjoint quarks in the large N limit in 
higher dimensions [0, Q where it has a substantially more complicated structure, although, 
even there, a solution of some special cases of the model are relevant to the deconfinement 
transition of three and four-dimensional Yang-Mills theory 0. 

In the present paper, we shall concentrate on solving a more general version of the one- 
dimensional model than has previously been considered and elaborating on the properties 
of the solution. The one- dimensional case has the advantage of being directly related to a 
continuum field theory, the heavy quark limit of 1+1-dimensional quantum chromodynamics 
(QCD). Part of the motivation for this work is to study the possibility of a confinement- 
deconfinement phase transition at high temperature or density in a field theory which has 
some of the features of QCD, i.e. with similar gauge symmetries and where interactions are 
mediated by non-Abelian gauge fields. QCD exhibits confinement at low temperature and 
density with elementary excitations being color neutral particles - mesons and baryons. On 
the other hand, at high temperature or high density it is very plausible that the dynamical 
degrees of freedom would be quarks and gluons - which would form a quark-gluon plasma, 
rather than the mesons and baryons of low temperature nuclear physics. At some intermediate 
temperature or density there should be a crossover between these two regimes. There are few 
explicitly solvable models where this behavior can be studied directly. Previous to the model 
of 0, Q , known explicit examples of phase transitions in Yang-Mills theory were not associated 
with confinement, but were either lattice artifacts unrelated to the continuum gauge theory 
or were associated with topological degrees of freedom in Yang-Mills theory on the sphere 
and cylinder p, Q and demark a finite range of coupling constants within which the gauge 
theory resembles a string theory. 

In finite temperature Yang-Mills theory (or QCD with only adjoint quarks), confinement is 
thought to be governed by the realization of a global symmetry which is related to the center of 
the gauge group and implemented by certain topologically non-trivial gauge transformations 
that appear only at finite temperature ||^, |10|, 0. The Polyakov loop operator is an order 
parameter for spontaneous breaking of this center symmetry and yields a mathematical way 
of distinguishing the confining and deconfined phases. When fundamental representation 
quarks are present, the center symmetry is broken explicitly and the Polyakov loop operator 
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is no longer a good order parameter for confinement. Whether, in this case, a mathematical 
distinction of confined and deconfined phases exists, and indeed whether there is a distinct 
phase transition at all, is an open question. 

Here, we will consider a toy model which resembles two-dimensional QCD with heavy 
adjoint and fundamental representation quarks. It could also be thought of as the heavy 
quark limit of dimensionally reduced higher dimensional QCD where the adjoint particles are 
the gluons of the compactified dimension which get a mass (similar to a Debye mass) from 
the dimensional reduction, and the fundamental representation particles are the quarks. We 
will solve this model explicitly in the large-A^ limit. The model with only adjoint quarks 
was solved in refs. 0, ^ and it was found that the explicit solution has a first order phase 
transition between confining and deconfining phases. These phases could be distinguished 
by the expectation value of the Polyakov loop operator, which provided an order parameter 
for confinement in that case. In this paper, we shall add fundamental representation quarks. 
Then, as in QCD, the center symmetry is explicitly broken and the Polyakov loop is always 
non-vanishing. We nevertheless find that the first order phase transition persists when the 
density of fundamental quarks is sufficiently small. When the density of fundamental quarks is 
increased until it is comparable to the density of adjoint quarks, the phase transition becomes 
a second order one. When the fundamental quark density is increased further, the phase 
transition is third order. 

It is tempting to identify this third order phase transition as the vestige of the deconfine- 
ment transition. Indeed, we shall present arguments for this. We shall find a mathematical 
way to distinguish the phases on either side of the third order transition based on the behavior 
of Polyakov loops with high winding numbers. We conjecture that the generalization of this 
argument to higher dimensional systems would provide a mathematical characterization of 
confinement in finite temperature QCD with fundamental representation quarks. 

Another motivation of the present paper, as well as Ref. P, |^ is to study a suggestion 
by Dalley and Klebanov and Kutasov [13| that l+l-dimensional adjoint QCD would be 
the simplest gauge theory model which exhibits some of the stringy features of a confining 
gauge theory. It is a long-standing conjecture that the confining phase of a gauge theory 



can be described by a string theory \T^\. There are only two cases where this relationship 
is well understood, two-dimensional Yang-Mills theory [|T^, [T^, |T^ and compact quantum 
electrodynamics |]TH]. At low temperatures l+l-dimensional adjoint QCD is confining in the 
conventional sense that quarks only appear in the spectrum in color neutral bound states. 
This is a result of the fact that, in one dimension, the gluon field has no propagating degrees 
of freedom and therefore it cannot form a color singlet bound state with an adjoint quark. As 
a result, the quantum states are color singlet bound states of two or more adjoint quarks. The 
spectrum contains an infinite number of families of multi-quark bound states which resemble 



asymptotically linear Regge trajectories [O, O, |T9| and, for large energies, the density of 



states increases exponentially with energy |2y]. This implies a Hagedorn transition |2T[ 
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at high temperature. Kutasov [|T3| supported this view by using an argument originally 
due to Polchinski that a deconfinement transition occurs when certain winding modes 
become tachyonic at high temperature. This behaviour was a feature of the explicit first 
order deconfinement transition found in the model considered in [0, 0]. That model, which 
coincides with the heavy quark limit of adjoint QCD, is effectively a statistical mechanical 
model for strings of electric flux, with quarks attached to their ends. 

The large- expansion of two-dimensional adjoint QCD has the same complexity as the 
large- iV expansion of a higher dimensional Yang- Mills theory and the leading order, infinite- 
limit cannot be found analytically In fact, the dimensional reduction of three or four- 
dimensional Yang Mills theory produces two-dimensional QCD with massless adjoint scalar 
quarks, so the combinatorics of planar diagrams is very similar. 



1 . 1 Overview 

This paper is organized as follows: In Section 2.1 we identify the gauged principal chiral 
model which corresponds to the gas with sources in various representations. In Section 2.2 the 
quantum mechanical formulation of this unitary matrix model is analyzed. This is followed 
by a section (2.3) where we rewrite the model in terms of collective variables (eigenvalue 
density), which are convenient for analyzing the large- iV limit. 

Parametrized solutions to the collective field equations are given in Section 3.1, and the 
parametrized free energy and its derivatives are obtained in 3.2. We also show that they 
give rise to a first order differential equation the free energy has to obey. In Section 3.3 we 
establish the existence of a third order line in the phase diagram, and compute the point where 
it terminates. Using numerical techniques we show in 3.4 that the critical line continues from 
that point first order line. 

In Section 4.1 we derive the contributions of the gauge field and the sources respectively 
to the energy density. This is followed by a section where we establish that higher windings 
of Polyakov loop operators behave differently in the two phases (4.2). This allows to discuss 
the phase diagram using energy densities and higher winding loop operators (4.3). The paper 
ends with a summary and outlook (4.4). Some of the more technical calculations are given in 
an appendix. 



2 Formalism 
2.1 Effective action 

The partition function of 1-1-1- dimensional Yang-Mills theory at temperature T and coupled 
to a number K of non-dynamical quarks at positions xi . . . , xk in representations . . . , Rk 
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of the gauge group is obtained by taking the thermal average of an ensemble of Polyakov loop 
operators 

Z[T; xi, . . . , xk; Ru...,Rk] = I dA, e"^!^! f[ TrPe^/o'"^-<' , (2.1) 

i=i 

where the Euclidean action is 

S[A] = fj^ dr j dx^i:i{F,,{r,x))\ (2.2) 

the gauge fields are Hermitean NxN matrix valued vector fields which have periodic boundary 
conditions in imaginary time 

A^{t,x) = A^{t + 1/T,x) , 

and the field strength is 

Ffiu = d^Ai, — d^A^ — i [A^, Ay\ . 

The gauge field can be expanded in basis elements of the Lie algebra A^^ = A^^T^. with T^. 
the generators in the representation For concreteness, we consider U{N) gauge theory 
and denote the generators in the fundamental representation as with a = 1, . . . , A^^. They 
obey 



rpa rph 



if'^'^T" , (2.3) 



normalized so that 

Tr T"T^ = -^-^^ (2.4) 



and with the sum rule 



a=l ^ 



We remark that group elements g^^ in the adjoint representation are related to the funda- 
mental representation matrices g by 

^gAdyb ^ 2 Tr(^^T"^T*) . (2.6) 

The expression ( |2.1| ) can be obtained by canonical quantization of l+l-dimensional Yang- 
Mills theory with Minkowski space action coupled to some non-dynamical sources 

S = - J dtdx ^TtF^^F^"" + source terms . (2.7) 
In the following, we will review an argument for representing the partition function of Yang- 



Mills theory as a gauged principal chiral model which was first given in and which was 



generalized to the case of Yang-Mills theory with sources in [0, ^ . As is usual in canonical 
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quantization of a gauge theory, the canonical conjugate E{x) of the spatial component of the 
gauge field (which we denote by A{x)), is proportional to the electric field, 



and obeys the commutation relation 



A\x),E\y)\ = zS'''S{x-y). (2.8) 
The Hamiltonian is 

H = jdx^-Y.{E''{x)f . (2.9) 

a=l 

This Hamiltonian must be supplemented by the Gauss' law constraint equation which is the 
equation of motion for Aq following from ( |2.7| ) and which contains the color charge densities 
of the sources 

g^{x) = (^E^ix) - r''^A''{x)E%x) + yn S{x - xA ~ . (2.10) 
\dx f^i ' J 

Here, the particles with color charges are located at positions xi, . . . ,xk- T^. are generators 
in the representation Ri operating on the color degrees of freedom of the i'th particle. 

There are two options for imposing this constraint. The first is to impose another gauge 
fixing condition such as 

A ~ 0, 

and to use the constraints to eliminate both E and A. The resulting Hamiltonian is 

H=Y. -^T%.®n^^^-^^ , (2.11) 

i<j,a 

which was considered in Ref. [|l|. It is the energy of an infinite range spin model where the 
spins take values in the Lie algebra of U{N). 

The other option, which makes the closest contact with string dynamics, is to impose the 
constraint ( 2.10 ) as a physical state condition, 

g^{x) ^phys = . 

To do this, it is most illuminating to work in the functional Schrodinger picture, where the 
states are functionals of the gauge field, ip[A] and the electric field is the functional derivative 
operator 

The time-independent functional Schrodinger equation is 

(2 £2 \ 

"y|5(M^ = £^"-[A;x„...,Xk] • 
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Gauss' law implies that the physical states, i.e. those which obey the gauge constraint ( p.lO| ), 
transform as 



where 

A3 = gAg^-igVgK 

is the gauge transform of A. 

For a fixed number of particles, the quantum mechanical problem is exactly solvable. For 
example, the wavefunction of a fundamental representation quark- ant iquark pair is 

^nA;xux,] = (Ve^^^''^^'^y , (2.12) 

where the path ordered phase operator represents a string of electric flux connecting the 
positions of the quark and anti-quark. The energy is ^-^\xi — X2I. For a pair of adjoint 
quarks, the wavefunction is 

^"^[A;xi,X2] = Tr fr^Pe'^-^' . (2.13) 



with energy ^^|x2 — Xi\. These energy states are identical to what would be obtained by 
diagonalizing the 'spin' operators in the gauge fixed Hamiltonian ( |2.11| ). 

Note that the wavefunctions ( p.l2|) and (p.l3|) are not normalizable by functional integra- 



tion over A. This is a result of the fact that the gauge freedom has not been entirely fixed, 
so that the normalization integral still contains the infinite factor of the volume of the group 
of static gauge transformations. 

In general, for a fixed distribution of quarks, a state-vector is constructed by connecting 
them with appropriate numbers of strings of electric flux so that the state is gauge invariant. 
The number of ways of doing this fixes the dimension of the quantum Hilbert space. If the flux 
strings overlap, the Hamiltonian can mix different configurations, so the energy eigenstates 
are superpositions of string configurations. However, this mixing is suppressed in the large- iV 
limit (i.e. the strings are non-interacting) and any string distribution is an eigenstate of the 
Hamiltonian with eigenvalue (e^A^/4) x (total length of all strings). 

We shall study the thermodynamics of this system by constructing the partition function. 
We work with the grand canonical ensemble and assume that the quarks obey Maxwell- 
Boltzman statistics. The partition function of a fixed number of quarks is constructed by 
taking the trace of the Gibbs density e~^/"^ over physical states. This can be implemented 
by considering set of all states in the representation of the commutator (|2.8| ), spanned by, for 
example, the eigenstates of A"-{x) and an appropriate basis for the quarks 
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Projection onto physical, gauge invariant states is done by gauge transforming the state at 
one side of the trace and then integrating over all gauge transformations (and then dividing 
by the infinite volume of the gauge group) PBl . The resulting partition function is 



Z[T/e2;xi,...,xx] = J[dA][dg] (A| e"^/^ Tr . . . Tr (^''"(a;^) , (2.14) 

where [dg] is the Haar measure on the space of mappings from the line to the group manifold 
and [dA] is a measure on the convex Euclidean space of gauge field configurations. The 
expression (|2.14| ) is identical to (|2.1|) with the Polyakov loop operator is the trace of the 



group element g{x) in the appropriate representation. 

In going over to the grand canonical ensemble the first step is to integrate over all particle 
positions. We then multiply by the fugacities for each type of charge: a factor of for 
each quark in representation R. To impose Maxwell-Boltzmann statistics, we divide by the 
factorial of the number of quarks in each representation. We then sum over all numbers 
of quarks in each representation. This exponentiates the fugacities, resulting in the grand 
partition function 

Z[T/e',XR] = J[dA][dg]e-'^-^'''^\ (2.15) 

where the effective action is 

exp{~S,f,[A,g]) = {A\e-''/^\A'^)exp(^Jdx A^Tr (?^(x) j , (2.16) 

and the summation in the exponent is over all the irreducible representations of U{N) we 
want to consider. The Hamiltonian is the Laplacian on the space of gauge fields. The heat 
kernel obeys the equation 

with the boundary condition 

Jim^(A|e-^/^|A^) = 5{A-A<^). 

These equations are easily solved by a Gaussian - divided by a T-dependent constant: 

{A\ e-^l^ \A^) ~ exp J dx ^ Tr {A - A^y^ . 

We see that the effective theory is the gauged principal chiral model with a potential 
energy term for the group- valued degrees of freedom, 

S,s[A,g] = jdx (^Ti\Wg{x)+i[A{x),g{x)]'\ -Y.{^R^^ 9\x)) \ . (2.17) 



R 
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Note that we have introduced the couphng constant 

7 - (2.18) 

When we analyze the hmit — >■ oo, we will tune such that 7 is constant. Moreover, we 
assume that the fugacities \r are scaled such that all terms in the action ( |2.17| ) are of order 
N\ 

The potential energy term in the effective action, 

V{g) ^ -Y.^nTi(g\x)) , (2.19) 

is the expansion of a local class function of the group element g{x) (one which obeys V{g) = 
V{hgh~^) for h G U{N)) in group characters with coefficients Xr. The characters 

XRig) = Tr((7^(x)) , 

form a complete set of orthonormal class functions of the group variable, with inner product 

[dg]x*Ri9)xB'ig) = 5r,r' ■ 



Here [dg] is not a functional integral measure, but is the Haar measure for integration on 
U{N). From the potential, we can find a fugacity by 



= -J [d9]x*R{9)V{9) ■ 
By tuning the fugacities appropriately, we could obtain any local invariant potential. 

The effective action ( |2.17| ) with all A^^ = was discussed by Grignani et.al. [^] and was 
solved explicitly in the limit ^ 00 by Zarembo [^]. The model with X^d 7^ (with adjoint 
quarks) was solved in Refs.pl and The effective action (|2.17| ) is gauge invariant. 



S,^[A,g] = S,^[A\hgh^] . 
It is also covariant under the global transformation 

S,s[A,zg,XR] = S,s[A,g,z-''^^'''>XR]. (2.20) 

where 2; is a constant element from the center of the gauge group, which for U{N) is ~ 
f/(l) and would be the discrete group for gauge group SU{N). Here, Ci(-R) is the hnear 
Casimir invariant of the representation R, which is the number of boxes in the Young tableau 
corresponding to R. When the gauge group is SU{N) and the only non-zero fugacities are 
for the zero 'N-ality' representations, i.e. those for which Ci{R) = mod N, there is a 
global Zjv symmetry. For gauge group U{N), this occurs only when all representations with 
non-zero fugacities have equal numbers of quarks and anti-quarks. The fugacities of other non- 
symmetric charges, can be thought of as an external field which breaks the center symmetry 
of the system explicitly. This situation is akin to the effect of an external magnetic field on 
a spin system. 
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2.2 Matrix quantum mechanics 



If we re-interpret x as Euclidean time, the partition function that we have derived has the form 
of a Euchdean space representation of the partition function for matrix quantum mechanics, 
where the free energy is identical to the ground state energy of the matrix quantum mechanics. 
We can study the latter model by mapping the problem to real time r by setting x = it and 
A —iA. The action in real time is then 



QM 



dr (^^,\g + ^[A,gf -V{g) 



We remark that this action must not be confused with the action (|2.2|). Sqm is the action 



for a 0+1-dimensional problem (quantum mechanics), while ( |2.2| ) is the action for Yang Mills 
theory in 1+1 dimensions. This remark also holds for the Hamiltonian below. In order to 
avoid confusion, we label the quantum mechanical quantities with the subscript QM. 

The canonical momentum conjugate to the group valued position variable g is the Lie 
algebra element 

N 



n = —iig^g + g^Ag- A 
7 ^ 



and the Hamiltonian is 

Hqm = ^Trn^ - ^Trn(^t^^ _ ^) + v{g) . (2.21) 
The gauge field A plays the role of a Lagrange multiplier which enforces the constraint 

gllg^ - n ~ , 

and the Hamiltonian reduces to 



QM 



7 

2N 



TrH^ + V{g) 



We can expand the canonical momentum as 

H = ^H^T'^ 

a 

Then, the components satisfy the Lie algebra 

H^ H^ 



-T-g^ . 



(2.22) 
(2.23) 
(2.24) 



It follows that in the Schrodinger picture the components of the canonical momentum are 
represented as 

d d 
og ^ dgik 
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Denoted in components the constraint reads 



G" = (^TrT"c/TV - ^5"^^ n'' ~ . 

The constraint has no operator ordering ambiguity. It generates the adjoint action of the 
symmetry group 

The constraint can be reahzed as a physical state condition 

In the representation where states are functions of this imphes that the physical states are 
class functions 

where h e U{N). This means that the physical states are functions of the eigenvalues of g. 
In a basis where g is diagonal, 
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the wavefunctions are functions of ccj, 

^phys("l, • • • , CUn) = V'phys(«l, • • • , ttj + 27r, Q;jv) 

Denoting the gauge group Laplacian in components 

JV2 



(2.25) 



(2.26) 



(2.27) 



o=l 



the Hamiltonian reads 



H, 



QM 
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A + V{g) . 



Since the potential V{g) is also a class function and depends only on the eigenvalues, when 
operating on the physical states, the Hamiltionian can be expressed in terms of eigenvalues 
and derivatives by eigenvalues 



Hqm — 



where 



J (a) = 2 sin -(ctj — aj) 



J{a) 



i<j 



(2^)^(^-^)/^^.4 



(N-l)/2 ' 
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and 

i<j 



The physical states must by symmetric functions of aj. (There is a residual gauge in- 
variance [^| under the Weyl group which permutes the eigenvalues and the physical state 
condition requires that the physical states be symmetric under these permutations.) The 
normalization integral for the wavefunction is 

[d9]^H9)^{9) = 1 • 

Since the integrand depends only on the eigenvalues of g, It is convenient to write the Haar 
measure as an integral over eigenvalues of g with the Jacobian which is the Vandermonde 
determinant, 

\lldai)\J{a)\'^tlj\a)^{a) = 1. 



The Hamiltonian and inner product have a particularly simple form when we redefine the 
wavefunction as 

ijj{ai, . . . , On) = J (a)V^(ai, . . . , on) ■ 

Since J is antisymmetric, ^ is a completely antisymmetric function of the eigenvalues, which 
we can think of as the coordinates of fermions. The Hamiltonian is that of an interacting 
Fermi gas 

(E - N{N' - 1)/12^ + i^ia) = Si,{a) . 



This correspondence of a c=l matrix model with a Fermi gas was first pointed out in Ref. |29 . 



2.3 Large N: Collective variables 

In this section we shall examine the collective field formulation of the large- iV limit of the 
theory that we discussed in the last subsection The Hamiltonian obtained in the last 

subsection reads 

Hqm = 7^E(n'^)' + Via), (2.28) 

a=l 



with (compare (|2.19| )) 



Vig) ^ -J2\nTv(g 



-Ml 



T 



R 

It was shown (compare ( |2.25| ),( P^ )) that the wavefunction depends only on the eigenvalues 
e*"-' of g and thus the density of eigenvalues 

1 ^ 

i=l 



11 



completely characterizes the properties of the system. Interpretation of the physics of the sys- 
tem at large N is more convenient when one considers the Fourier transform of the eigenvalue 
distribution 

p{0,r) = - + ^Ecn(r)e-^ (2.29) 



where we have defined 
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c„(r) = -Tr(7"(r) , c_„(r) = c„(r) . 

We now turn our attention to developing the collective field theory formulation of the (thermo-) 
dynamical problem given by the Hamiltonian ( p. 2^ ). Since the wavefunction depends only on 
the eigenvalues of g, we would like a Hamiltonian equivalent to ( p.28| ) but written in terms of 
the eigenvalue density p and a conjugate momentum 11. At large we will find this Hamil- 
tonian and write equations of motion for p and H. So far we have not imposed any restriction 
on the potential V{g), but from now on we assume, that it can be expressed as a functional 
of the eigenvalue density p{0). In particular the potential we are going to analyze below will 
have this property. 

The canonical momentum operates on the wavefunction as 

= A? /^^Ee-^'''^Tr(T"^^)-47jr 

and the Laplacian ( p.27|) is 

/ 1 /• / ^ \ S 



8%^^ J ^ ^ Sp 9)Sp e') ^^^^ 



which can be written as 



= -T^r I d9p{9) \ ( 4t47IvV - N'V f d9'p{9') cot ( ^—^] ^^-^ \ ^[p] 




06 6p{e) 



where 



v{e) = —V j de'p{e') cot 

V indicates principal value integral. 
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The transformation of the wavefunction 



^[p] = vi[p]exp(^-^ ldede'\nsm^-^-^p{e)p{e')^ , 

transforms the derivative in the Schrodinger equation so that it has the form 



The second term on the left-hand-side of this equation has a simple form. In Appendix A.l 
it is shown that it gives rise to a term which is cubic in the density. Thus, the Schrodinger 
equation has the form 

ihh (mmJ- T^'w) + ''w} *i = ^ *i . (2.31) 

up to an overall constant. 

The large- limit is dominated by the eikonal approximation. In this approximation, we 
make the ansatz for the wavefunction 

i^[p] = exp {iN'S[p]) . 

The eikonal, S then obeys the equation 



m 8 



hri^m) - T^=(«) h ^^'^w ^ ^ . (2.32) 



Here, we have ignored a term which is of subleading order in A^^. We have also assumed that 
V[p] will be of order A^^ (compare Section 2.1) and that the natural magnitude of the energy 
eigenvalue is of order A^^. 

To solve this equation for the ground state, we must find its minimum by varying p and 
the canonical momentum 

n = 5S/5p, 

subject to the condition that p is normalized. This leads to the equations of collective field 
theory 



dr' ' ' SU{T,e) ' 

-y{r,e) = 



9t ^ ' ^ d0 Sp{T,e) ' 

where 

v{T,e) ^ ^n(r,^). 
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Taking the derivative of the second equation with respect to 9 ehminates a Lagrange multipher 
which must be introduced on order to enforce the normahzation condition for p. Using ( |2.32| ) 
one finds 

d 79,. 

|„,2|(„.,.V).^||;.W ^ 0. (2.33, 

It is interesting to note that these are nothing but Euler's equations for a fluid with equation 
of state P = 7r^p^/3 on a cyhnder with coordinates {6, r). The inclusion of a potential V{6, r) 
corresponding to non-Abelian charges is equivalent to subjecting the fluid to an external force 
which is derived from V{6, r). 

We shall use these equations in the next section where we analyze the large- limit of a 
mixed gas of adjoint and fundamental charges. 



3 Free energy and critical behaviour 



3.1 Static solutions to the collective field equations 



In this section we will find static solutions to the collective field equations (p.33|) . The most 
simple potentials involve only the lowest representations, the fundamental, its conjugate and 
the adjoint. We shall consider a slight generalization of these and use powers of the lowest 
representations to include multiple windings of the Polyakov loop operator. Our potential 
readsQ 

V{g) ^ -£(K„iVTr(^7«) + /s:„iVTr(((7t)«)+A„|Tr^7"|2) , (3.1) 



n=l 



where we made use of (compare ( p.5|) , ( \2-Q[ 

Tr (g^^x)) 



|Tr g'^ix) 



to relate the trace in the adjoint representation to the trace in the fundamental representation. 
The couplings for the fundamental representation charges (and their conjugates) were chosen 
to scale ~ A^, to make the potential of order N"^. 

The potential ( |3.1| ) indeed can be expressed as a functional of the eigenvalue density ( |2.29|) . 
The collective field Hamiltonian ( |2.32|) then reads 



QM 



^ I de 



p{0) He)f + \p\o) 



(3.2) 



n=l 



de p{e)e 



ind 



+ Kn de p{e) e*"' + Kn de p{e) e 



-inO 



1_ 

96 



^It should be remarked that parts of this section can easily be extended to more general potentials. 
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In order to maintain correspondence with the original version of the Hamiltonian ( |2.28| ) we 
subtract the constant 7/96. It sets the energy scale such that the free energy vanishes in the 
confined phase of the model with only adjoint charges (see below). 

The corresponding collective field equations ( p.33|) read 



dv 7 (9f ^ 7r^7 dp^ d 



where c„ are the x-dependent Fourier coefficients of p as introduced in ( p.29| ). Note that we 



also performed the change of variables, r — > —ix and f — > in these equations in order 
to invert the Wick rotation performed at the beginning of Section 2.2 prior to canonical 
quantization. 

We will only consider real, static solutions of the non-linear equations (p.4|) , that is, where 
p{0,x) = Pq{9) and the velocity v vanishes identically. Consequently 



^Q\ j 2 J-^V-E + E(A„c_„ + K„)e*"^ + E(AnC„ + Kn)e where p is real , . 
P'^^^ = \ otherwise " ^'-'^ 

The constant of integration E has physical interpretation as the Fermi energy of a collection 



of fermions |^9[| in the potential V[p] and is fixed by the normalization condition 

1 = Jde po{e) . (3.6) 

Here it is more convenient to express the c„ in terms of po (compare (|2.29|) ) 

d9 po(^)e'"^ . (3.7) 



The real support of the function po{0) is the positive support oi A = E + J2{^nC-n + i^n)^'^"'^ + 
J2{^nCn + ft;„)e~*"^. The zeros of A define the edges of the eigenvalue distribution and when 
these zeros condense one has critical behaviour in the observables of the model as in general 
Hermitean and unitary matrix models. 

3.2 A differential equation for the free energy 

In this subsection we compute all first order derivatives of the free energy and show that they 
obey a differential equation of the Clairaut type. 

Inserting the static solution ( ^.51 ) in ( |3.2| ) we obtain for the free energy 



1 11^ 7 

^(^QA/) = f = - - ^ [A„|c„P + 2( K^c^ + K„c_„ )] - — . (3.8) 
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Note that / is the leading coefficient (0(A^^)) of the energy for the matrix quantum mechanics 
problem, but for the quark gas problem plays the role of the leading coefficient of the energy 
density. 

Deriving the expression ( |3.8|) with respect to Aj for some fixed J and using derivatives of 
Equations ( p.6|) and ( |3.?| ) with respect to the same parameter one obtains 

df 



dX:, 



-|cj| 



(3.9) 



We use the notation d/d\j to indicate that also the c„ and E which implicitly depend on Aj 
are derived with respect to this coupling. Similarly one can show 

df 



and 



d'j 



dK.j 



1 1 oo 

37 37 L 



-cj 



1 

96 



(3.10) 



(3.11) 



It is interesting to notice, that combining Equations ( |3.9| ) - ( |3.11| ) gives rise to a ffist order 
differential equation of the Clairaut type 



7- 



df 



+ 

n=l 



df 



df 



+ Kr 



df 



d'y ~, L dXr, dKr, dUr 

This differential equation has general solutions of the form 



f 



(3.12) 



/(7, A„, Kn) 



7 7 



where F is some arbitrary smooth function. This result shows that the parameter 7 is not 
driving the physical properties of the model, but rather sets the energy scale. The differential 
equation ( |3.12|) gives no further restrictions on the function F and another analysis will be 
adopted in the next section. However, when finding the physical interpretation of the phase 
diagram in Section 4.3, the differential equation (|3.8|) will be a valuable tool. 



3.3 Regime of the third order phase transition 

We now restrict ourselves to the case of only one pair of non-vanishing couplings Aj, kj 7^ 0. 
Furthermore it is sufficient to consider kj real, since an eventual phase of nj can always be 
removed by using the covariance ( f2.2(J| ) of the action and Haar measure [dg] under transfor- 
mations by a constant element of U{1). 



In the form of (3^) it is evident we need to solve simultaneously for the normalization 
condition ( |3.6| ) and the Fourier moment ( p.7|) in order to have a self-consistent solution of the 
saddle-point equations. We begin by introducing an auxiliary complex parameter 

rj e'"'-' = Aj c_j + Kj , (3.13) 
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and rescaling the Fermi energy as 

E = 2firj. 

With this notation the normahzation and moment equations are respectively 




cje 



ljf3 



2 M(/i) , 



(3.14) 



(3.15) 



where we have defined the integrals 

2 



M(/i) 



+ cos{J6) H(/i + cos( J6')) de , 
- r cos J9JfL + cos{J9) H(/i + cos( J^)) d9 . 

TT J — n 



(3.16) 



Here H(..) denotes the step function. A simple transformation of the integration variable 
shows that /(/i) and M(/i) are independent of J. Thus J enters only as the subscript of the 
parameters. For notational convenience we abbreviate 



A., 



A 



Kje 



K 



(3.17) 



We remark that /(/i), M{fi), and thus cje^'^^-' are real. Eliminating the moment cj from (|3.15|) 
by using the definition (|3.13| ) we obtain 



K 

7 



- - -/(/.)M(/.) 

4 7 



(3.18) 



This family of lines in the k, A-plane parametrized by yU represent a necessary condition which a 
solution of the normalization and moment equations ( p.l5| ) must obey. From the last equation 
it is obvious that also the product k = Kje~^'^^ is real (thus e*"^^ is just a sign). It occurs as 
a natural parameter when rewriting the free energy in terms of /(/i) and M(/i) (use ( p.8| )) 



/ 



7 



12/(/.)= 



2/1 



M(/i) 



1_ 
96 



(3.19) 



where we have eliminated A using the necessary condition (|3.18|) . Also the first derivative of 
the free energy with respect to k can be expressed conveniently in terms of /(/i) and M(/i) 



-2 



(3.20) 



Remember that we restricted ourselves to kj real, and thus we encounter a factor 2 compared 
to 3.101) , since a real kj is the same for both terms cj and c_j in the potential ( p.8|) . 

With the parametric solution ( ^.8[ ) at hand we turn our attention to establishing the 
critical behaviour in this model. In the appendix it is shown that the first derivatives of 
and M(/i) have non-analytic behaviour at /i = 1, hence the expression ( p.20|) suggests that 
the vicinity of /i = 1 is a natural place to look for non-analytic behaviour in the free energy of 
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our model. Using the explicit results ( [A.5| ) for /(I), M(l) and ( |3.18| ) we obtain the necessary 
condition for the critical (/i = 1) values of A and /t 



7 



TT 



512 



1 A-^ 



(3.21) 



Having identified a line in the phase space where we expect critical behaviour we will now 
proceed to establish the details of this critical behaviour. Following we begin by expanding 
about fi = 1 and the line ( |3.21| ) 



11 = 1 + e,e > 



1(1 + e)=h + 6I 



M(l + e) = Mc + 



(3.22) 



(/c = -^(1), Mc = M{1)) and analyze the variation of k around k'^ while keeping A fixed 



A = A^ 



(3.23) 



Due to the nontrivial support of the integrands in ( p.l6|) in principle one has to distinguish 
the cases /x > 1 and /i < 1; (see the appendix for details). In order to keep the formulas 
simple, we explicitly analyze only the case /i > 1 as given by (|3.22| ), (p.23|) . The case /i < 1 
can be treated along the same lines and we denote the corresponding results in the end. 

The expansion now consists of two steps. We first expand the necessary condition ( |3.18| ) 
at /i = 1 to obtain the relation between the variation 5k and e. In the second step we expand 
the right hand side of ( ^.2(j| ) at = 1 and use the result of step one to express the variation 
of df/dK in terms of 5k. The latter result can then be used to analyze eventual singular 
behaviour of higher derivatives of the free energy. 

Expanding the necessary condition (|3.18|) and using ( 3.21 ) we obtain for the variation of 
K to lowest order 



5k 

7 



1 51 



2(/c 



7 Ic 



5M 



51 

T 



-7r2 4 A^ 
256 ^ 3" 



A^£ 

72 ' 



(3.24) 



where in the last step we made use of the relation ( [A.9|) between the variations 51 and 5M 
and inserted the explicit results ( |A.5| ) for 1^ = /(I) and Mc = M{1). Using the result ( |A.8 
for 51 /Ic to lowest order we obtain 



5k 



where we introduced the abbreviation 



a 



— e\n{e) a 



-fii^ 4A" 
256 ^ ~3~ 



Inverting equation (|3.24|) (again taking into account only the leading order) gives 

— a~^5K In (a~^5K^^ 



(3.25) 



(3.26) 



(3.27) 
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This equation is the relation between the variation 6k and e which is imphed by the necessary 
condition ( p.l8| ). In the final step we expand the derivative of the free energy ( p.20| ) at = 1 
and use the result p.27| ) to obtain the variation of the derivative in terms of 6k. Expanding 



(lOOl ) gives 



dK 



Me 

-2 — 



1 + 



6M 61 



Mc 



2 
3 



(3.28) 



Using ( 3.27| ) we obtain 



dl 
dK 



-- + -a ^6k + cr ^6k In f 0" ""^^kV 
3 3 I \ ) 



(3.29) 



We remark, that the case /i < 1 with expansion fi = 1 — e, e > changes only the sign of the 
argument of the logarithm. Differentiating the last result with respect to 6k establishes the 
singular behaviour of the third derivative of the free energy with respect to k. Thus we find 
a third order phase transition for /i = 1. The critical line is a straight line given by ( |3.21| ). 

It is important to notice, that at (see Equation (|3.24|) ) 



7 



37r^ 
1024 



(3.30) 



the leading term in the expression for 6k vanishes. Equation ( p.24|) is reduced to the simpler 
relation 



5k 

7 

At this point the expansion of df /dK gives 



—e 



27 



dK 



2 
3 



-e ln(£) 



2 2^/2^ 
-- + — — OKm ——6k 

3 3A^ V 



2 ^ 
+ ^6k 



(3.31) 



(3.32) 



Again the case /i < 1 differs only by the sign of the argument of the logarithm. Differentiation 
with respect to 6k shows, that the phase transition has turned to second order at that point. 
Using ( p.21| ) one can compute also the k/j coordinate of the second order point giving k/j = 
71^/1024, A/7 = 37r^/1024. In fact the more global analysis in the next section will show, that 
the third order line terminates at the second order point /t/7 = 7r^/1024, A/7 = 377^/1024, 
and continues as a first order line. 



3.4 Regime of the first order phase transition 

As pointed out at the end of the last section at the point k/'j = 7r^/1024, A/7 = 37r^/1024, 
the third order transition along the ^ = 1 line ( p.21| ) changes to second order. This unusual 
behaviour requires further investigation which we will carry out in this section. To begin, 
a graphical analysis of the phase diagram is most useful and in Figure we plot a number 
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Figure 1: Plot of the lines (3.18) for n ranging from 0.4 (upper right corner) to 75 (line at 
the extreme left). The region of overlapping lines corresponds to a region of first order phase 
transition. 



of representatives of the family of lines (|3.18| ) for a range of values of /i. It is clear that in 
most of the n, A-plane points are in a one-to-one correspondence with values of the parameter 
/i. This correspondence breaks down though in a small region near the A/7 axis between 
A/7 = 0.05 and A/7 = 0.0625. Due to the behaviour of the slope and intercept in the linear 
equation (|3.21| ) lines begin to overlap for increasing fi starting at /i ~ 1 and continuing as 
H — >■ 00. In this overlap region the phase diagram is folded at the vertex n/'j = 7r^/1024, 
A/7 = 37r^/1024, and each point falls on three different lines of constant /x. Consequently the 
system simultaneously admits three configurations with different free energies in this region of 
the phase space. This circumstance allows for a first order phase transition to develop along 
a line where the free energies of the different phases are equal. 

The edges of the triangular first order region in Figure |^ is given by a caustic of lines 
from the one parameter family ( p.l8| ) . The boundary is defined by the curve where the family 
of curves is stationary with respect to fi. This condition can be used with ( p.l8| ) to give a 
definition of the boundary caustic. The stationary condition can be solved with the parametric 
result 

K ^ 1 I'{f,)M{f,) + M'{fi)I{f,) 

7 4J2(^) M'(/x)/(/i) - /'(/x)M(^) • ^ • ^ 

As can be seen, the curve given by ( p.33| ) intersects the A/7 axis at two points: 0.057024 (yU = 
0.95324) and 1/16 (yU = 00) and reaches a singular maximum in the K/7 direction for /i = 1 
at the point K/7 = 7r^/1024. The end of this region of first order transitions agrees with the 
position of the second order transition point which was determined by the analysis of critical 
behaviour in the previous section. 
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Figure 2: Plot of the boundary of the multiple phase region. The boundary (solid curve) is 
given by a caustic of lines in the one-parameter family (3.18). The dotted curve shows the 
numerically determined first order line. 



Once one has determined the region where different phases can co-exist the next issue to 
address is that of the position of the hne of first order phase transitions where different phases 
have the same free energy. This hne can be determined for given K/7 by the simultaneous 
solution for the parameters /ii and ^2 of the pair of equations 

4 /(/ii)/(/i2) 7^ 

and 

Q[p{fx^) P{iX2)j 12U'M I'M) 7 /(/ii)/(/^2) ' ^' ^ 

Unfortunately, these equations are analytically intractable. Again we turn to a graphical 
analysis to gain further insight. In Figure |^ we plot the free energy of the system as a function 
of A/7 for different values of fixed K/7. From here it is easy to see a number of features of the 
region of first order transitions. Increasing /i traverses these curves in a clock-wise rotation 
so that free energy increases for small values of /x, intersecting the nearly horizontal large 
free energy. This intersection point is a graphical demonstration of the first order transition 
which occurs here as the model jumps from weak (/i < 1) to strong coupling (yU large). Each 
phase continues to exist after the transition point and may be reached by an adiabatic process 
until ending in cusps which mark the boundaries of the first order region in the A/7 axis. It is 
interesting to note that there is an energetically unfeasible intermediate "medium coupling" 
phase which connects the weak and strong phases. Hence, for fixed K/7 there exist three 
distinct configurations of the system for given A/7 in the region of first order transitions. 



(3.34) 
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Figure 3: Free energy j as a function of A/7 in the region of first order phase 
transitions. Each curve is plotted for fixed k/j which from right to left is k/'j = 
0, 0.0005, 0.001, 0.0015, 0.002, 0.0025 

4 Discussion 

We show that the whole phase diagram can be understood in terms of the energy densities 
of the fundamental charges the adjoint charges and the flux lines. Higher winding Polyakov 
loop operators will be used to characterize the phases of the model. 

4.1 Energy densities 

In this section we identify the expressions which should be interpreted as the energy densities 
of the fundamental charges the adjoint charges and the flux lines. 

Using the deflnition ( |3.8| ) for the leading (0(A^^)) coefficient / of the density of the free 
energy we can formally denote the partition function Z[T/e^, A„, defined in ( |2.15| ), ( |2.16D 
as (note that 7 is related to T/e^ via ( f^.l8D ) 

Z[T/e^,X^,^r.] ~ exp (-iV^L /(7)) . (4.1) 

L denotes the space-like extension which is infinite. Thus ( [4.1|) as it stands can only be 
understood formally, but could be made an exact identity by putting the system in a finite 
box. However, we show that derivatives of / have an interpretation as energy densities also 
in the infinite system. 
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Using ( p.l5| ), ( p.l6| ) we obtain 

d\nZ 



dT- 



df 



dT- 



where in the last step we made use of ([4.1|) . Using the definition for 7 (|2.18|) we obtain 



dl 

d'j 



1 (Hpiux) 

iV2 L ■ 



/ 2^^ _ 

e — = pFh 



(4.2) 



Pfiux is the leading coefficient (0(iV^)) of the contribution of the electric fiux to the energy 
density, measured in units of e^N/2. Note that is a proper unit for energy densities, 

since it is invariant in the limit N ^ 00 (compare ( p.l8|) ). The inclusion of the factor 1/2 is 
natural, since the quadratic Casimir operator behaves (for U{N)) as N/2 at large N. 

In the same way we can proceed to identify the contribution of the adjoint particles to the 
energy density^ 

Aj df 



N 



1 (HAd) I 2_ 

iV2 L / 2 



= pAd 



7 c?Aj 

Finally for the contribution of the fundamental charges and their conjugates we obtain 



7 dKj iV2 L ■ 
The differential equation ( p.l2|) gives 



/ 2^^ _ 

e IT = Pf.f 



(4.3) 



(4.4) 



PFIux — PAd ~ Pf,F 



7 



1 F 



,N 



N^L'^ 2 



(4.5) 



where F denotes the free energy of the system. Thus the density of the free energy is the 
difference of the contribution of the flux and the contribution of the charges. 

When interpreting the phase diagram it is convenient to express the densities in terms of 
ft/7, A/7, /(/i) and M(/x). Using (p]), ( glOD , ( CT) and we obtain 



1 



Pfi^ 



p + 



M{p) 



pAd 



Pff 



6/(/i)2 

A /M(/i)Y 

7 I Hp) J 

o^M(/i) 

7 I{p) ■ 



Hp) 



1 

96 



(4.6) 

(4.7) 
(4.8) 



The fact that ppiux does not depend on k./'J and A/7 is most remarkable. It proves that the 
lines of constant /i are the lines of constant flux. In the last section we showed that the lines 

^Again we restrict ourselves to only one (real) pair Aj, kj 7^ 0, although some of the results of this section 
can be immediately taken over to the case of several nonvanishing couplings. 
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of constant /i are responsible for the structure of the phase diagram. This imphes that the 
whole phase diagram can be understood by the abundance of flux. 



We remark that the necessary condition (|3.18|) gives rise to another quantity which is 



constant on lines with fixed /i. It is the linear combination 

4.2 Higher windings of Polyakov loop operators 

In addition to analyzing the phase diagram in terms of contributions of the electric flux and 
adjoint and fundamental charges to the free energy we also investigate higher windings of the 
Polyakov loop operators. As we will now demonstrate, their behaviour changes fundamentally 
at the line fi = 1. It can be used to characterize the difference between the regions /i > 1 and 
fi < 1 and the physics in these different phases. 

In the model with only the pair kj, Aj 7^ we will be concerned with the quantities 

PM = lim 1 < Tig''-' > = cuj. (4.10) 



Using (13), (|J|) and {^J§ we obtain 



cue'^" = ^ . (4,11) 



where we have defined 



Mkifi) = - r cos kO^fi + cos(e) H(/i + cos(e)) d9 . (4.12) 

Note that as in the cases of /(/i) and M[ii) (see Appendix A. 2), we have performed a transfor- 
mation of the integration variable which removes the dependence on J . The most important 
thing to note about this last equation is that these higher windings depend only on the pa- 
rameter /i and hence give a unified picture of different regions of the phase space as does the 
electric flux density ppiux- 

In Appendix A. 3 we show that the Mfc(/i) and thus the higher windings of the Polyakov 
loop operators Pfc(/i) undergo a drastic change of character at the point /i = 1. We find, that 
for /i > 1 they are exponentially suppressed with increasing fc, while for < ;U < 1 there is 
only power law decay. Our results are (use ( |A.15| ), ( |A.13| )) 



|Pfc(/i)| < (r(/i))'^v'2/i + r(/i)+r(/i)-i for > 1 , (4.13) 



with r(yu) = (1 + yU — ^/ jjfl — l)/2 < 1 for /i > 1. For /i = 1 we obtain 



= ip^' (4-14) 



24 



and for — 1 < < 1 it is shown in Appendix 3.1 that exponential suppression is excluded. 

In terms of a Fourier decomposition of the eigenvalue distribution po{9), this result is a 
comment on the smoothness of the distribution. For /i > 1, po{6) has support on the entire 
circle and is infinitely differentiable, whereas for n < 1 the zeros of po{G) with their infinite 
slope contribute higher frequencies to the Fourier series. 

Physically one would like to be able to associate the change in behaviour of the higher 
winding Polyakov loops to a change in character of the gas of fundamental and adjoint charges. 
As seen above, and as will be expanded on in the next section, the phase diagram may be 
broadly divided into two regions. We have proven that p = 1 denotes a line of third order 
phase transitions up to a point where a first order line develops and there is a competition 
between phases with fi < 1 and /i > 1. Thus the behaviour of the Pk{p) is different on both 
sides of the third order section of the critical line, as well as on the two sides of the first order 
section. 

It is interesting to notice, that the large-fc behaviour of the higher winding Polyakov loop 
operators Pk{p) can also be interpreted as the long distance behaviour of a gauge invariant 
two-point function. In this picture Pk{p) then corresponds to a Polyakov loop winding k- 
times around compactified time of the Euclidean model (compare Thus the Pk{p) can 

be considered as the expectation value of a gauge invariant two-point function with time-like 



separation fc/T, k = 1,2,3, .... For other models it was already proposed earlier to use 
the changing decay properties of such operators as a confinement criterion. 

It remains an interesting open question if the change in behaviour of the higher winding 
Polyakov loops persists in higher dimensions and for finite N. If it does, it would provide a 
powerful tool to study deconfinement transitions, in particular in lattice simulations, where 
the Pk{p) are much simpler to analyze than other criteria, as e.g. the area law of the Wilson 
loop. 

4.3 Physical interpretation of the phase diagram 

Figure ^ gives a schematic plot of the phase diagram. In the following we explore the physical 
behaviour of the system in various regions of the phase diagram. 

The line k = (purely adjoint model) 

In the case of only adjoint charges (all kn = 0) the action is invariant under transformations 
in the center of the gauge group, g ^ zg (compare ( |2.20| )). When this symmetry is faithfully 
represented we expect 

c„(x) = lim ^(Tr ^"(x)) = Vn ^ , (4.15) 



smce 



Tr ^"(x) ^" Tr ^"(x) . 
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Figure 4: Schematic picture of the phase diagram. The doted curve marks the first order part 
of the critical line. The solid curves above and below it are the boundaries of the area with 
two possible phases. They join at a point which shows second order behaviour. For larger n/'j 
we find a third order line (fi = 1) marked by a solid line. The numbers label four extremal 
corners of the phase diagram, which we discuss below. We also indicate the range of the 
auxiliary parameter /i. 



Prom ( p.3|) , ( p.4|) it is obvious that (|4.15|) which corresponds to constant eigenvalue density 
p = l/27r is a solution. The stability analysis performed in shows that this solution 
("confined phase") is stable for A/7 G [0,1/16]. Since the upper boundary of the multiple 
phase region is at A/7 = 1/16 we confirm this result. The confined phase at k/j = 
corresponds to /i = cxd (compare Fig. p. Inserting the results ( |A.10| ) for the large-/i expansion 
in ( |4.6| ), ( |4.7| ) and (|3.19| ) respectively, we find energy densities Pfiux = 0, pAd = and / = in 
the confined phase of the purely adjoint model. Also the higher winding Polyakov loops Pk{p) 
vanish as /i — 00. All this supports the picture of adjoint charges being bound pair-wise in 
singlets under center transformations. The phenomenon of adjoint charges bound in such 
"hadrons" is called confinement. 

For higher A/7 the center symmetry is spontaneously broken. This corresponds to a 
solution, where the coefficients c„ are non- vanishing as well as the densities Pfiux, PAd and the 
energy density /. Typical values (taken at A/7 = 0.07) are, ppiux = 0.017714, p^d = 0.020385, 
//7 = —0.002671. This "deconfining phase" has a nonvanishing overlap with the confining 
phase and gives rise to a phase transition which is of first order. The corresponding critical 
coupling can be computed numerically giving A/7 = 0.059250. This result is in agreement 
with the value given in 
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The line A = (purely fundamental model) 

The center symmetry of the action is exphcitly broken when there are fundamental charges 
(some Kn 0). In particular this is the case when we consider the model with only funda- 
mental charges (and their conjugate). Thus we cannot expect to find a confining solution 
with vanishing c„ which would correspond to vanishing densities Pfiux, Pf,f- 

For the purely fundamental model we have established the existence of a third order phase 
transition ai k/j = 7r^/512 (set = 1, A = in ( |3.18|) ). The phase with k/j smaller than the 
critical value is characterized by low densities Pfiux, Pf,f exponential suppression of the 
higher winding Polyakov loop operators Pk{p), while for k/'j greater than the critical value 
we find higher densities and power-like suppression of the Pk{p)- 

Using the expansion ( |A.10| ) for large fi, one can compute the behaviour of the densities at 
small k/'j. We obtain (expand ( p.l8| ) for large /x, to find the relation between k/'j and /i) 

PFiux = 16 (^^) + 0{iK/^Y) , p^,^ = 32 (^^^ + 0{{K/^y) . 

It is interesting to notice, that the two densities vanish with fixed ratio 2. This could be 
interpreted as a dilute gas of "hadrons" formed of pairs of fundamental charges and their 
conjugates bound together by flux. 



The third order line 

In Section 3.3 we established that at /i = 1 there is a third order phase transition. For A/7 < 
37rVl024 the critical line is given by Equation ( p.21| ). The third order behaviour terminates 



at k/'j = 7r^/1024, A/7 = 37r^/1024 becoming second order at that point, continuing as a line 
of first order transition. Using the values ( [A.5|) we obtain for the energy densities of adjoint 



and fundamental charges 

_ 2k _ 1A 

Pf,F — o~ ' PAd, — -- ■ 
37 97 

The energy density of the flux is constant along the third order line and given by (evaluating 
(13) at /i = 1) 

pFiux = 0.006718. 



The invariant linear combination ( [4. 91 ) assumes the value 

Pf^f + 2pAd = 0.012851 , 

on the critical line. 

The two phases above and below the third order line differ by the values of the energy 
densities, with higher values of density for large k/'j and A/7. In particular the fiux-density 
is higher for /i < 1, and lower for /i > 1. Thus we call the phase with p > 1 the "low density 
phase", while the other phase will be referred to as "high density phase". Since the energy 



27 



density of the flux ppiux vanishes entirely only for the purely adjoint model discussed in the 
last paragraph, we reserve the term "confined/deconfining phase" for this case. 

The most dramatic change is seen in the behaviour of the higher winding Polyakov loop op- 
erators We find exponential suppression with increasing k for /i > 1 which corresponds 
to the low density phase (small n/'j, X/'j). For < 1 (high density phase) we showed that the 
Pfc(/x) decay only power-like with k. On the critical line (/x = 1) we find \Pk{f^)\ = l/(4fc^ — 1). 
Together with the behaviour of the flux density ppiux, this allows to understand the whole 
phase diagram by the abundance of flux. 

The first order line 

The first order line comes about through the coexistence of two phases with equal free energy 
(compare Figure |^). The region where the phases coexist is shown in Figure 0. The two 
phases can be distinguished by the different contributions of flux and sources to the energy 
density (high density and low density phase). Typical values are (taken at K/7 = 0.002, A/7 = 
0.05203) 

Pfiux = 0.003328 , pAd = 0.002811 , pp^p = 0.000929 , 

for the low density phase which connects to the confined phase of the purely adjoint model 
(k/7 = 0,A/7G [0,1/16]) and 

Pfiux = 0.008999 , pAd = 0.007856 , Pf,f = 0.001554 , 

for the high density phase which connects to the deconfining phase of the purely adjoint model 
(fi;/7 = 0, A/7 > 0.057024). Although unlike for the purely adjoint case, for k/'J > we have 
no immediate interpretation using symmetries, since for K/7 > the center symmetry is 
always broken explicitly. However, the physical behaviour in the first order region is still 
reminiscent of the purely adjoint case. The values of all three densities are considerably 
enhanced in the high density phase. Since the numbers we quoted are taken for a point on 
the critical line, they give rise to the same free energy (Note that the free energy is not the 
sum of the three contributions, but rather related to them by ( |4.5| ).). 

It is important to notice, that also for the first order section of the critical line the higher 
winding Polyakov loops provide a criterion to distinguish the two phases. At the first order 
line a phase with p > 1 and a phase with p < 1 coexist. Thus also for the first order transition 
we find the characteristic exponential decay in the low density phase, but power-law behaviour 
in the high density phase. 

The second order point 

In the end of Section 3.3 it was established that the point k/j = 7r^/1024, A/7 = 37r^/1024, 
which separates the third order part of the critical line from the first order part, is of second 
order. At this point the energy density for the adjoint and fundamental charges respectively 
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are given by (insert /i = 1 and the coordinates of the second order point in ([4.7|), ([4 .81)) 

Tf_ _ Tf_ Pf,F _ ^ 

3072 153d PAd 

If one decreases k/7 from its value at the second order point, the adjoint charges start to 
dominate, and the ratio ppp/ pAd decreases - the first order region emerges. Conversely, with 
increasing /t/7 the fundamental charges start to dominate the system and Pf,f/ pAd increases. 

The existence of the phase transitions which we have discussed in this section is a bit 
surprising, considering the low dimensionality of the system. In fact, if we first consider the 
model with purely adjoint quarks as discussed in f^, there is a discrete symmetry under 
transformations in the center of the gauge group and the system is one- dimensional. Such a 
system cannot have a phase transition if N is finite. When is infinite, the center symmetry 
becomes U{1) and the group index behaves like another dimension, so the effective dimen- 
sionality of the model is two. In this case, the symmetry can be spontaneously broken (and 
there can be deconfinement) because the interaction potential has infinite range in eigenvalue 
space. This is also seen by considering the translation mode in eigenvalue space, which is 
the derivative by a of the eigenvalue density p. When the support of the density is a subset 
of the interval [— vr, vr] the derivative of p is not a normalizable wavefunction (divergence in 
the norm coming from the edge of the distribution) for the zero mode and thus it is not 
effective in restoring the symmetry. Similarly, we expect that, when we turn on the density 
of fundamental representation quarks, the first or third order phase transition happens only 
at infinite N . 

One can think of the infinite N limit here as something analogous to the infinite volume 
limit in the case of a ferro-magnetic system. In the latter, a phase transition in the math- 
ematical sense only occurs when the volume is infinite. However, in the physically relevant 
case where the volume is large but finite, the phase transition is clearly seen. Similarly, in the 
present case of large A^, the phase transition only occurs in the strict mathematical sense in 
the large limit. It should nevertheless be visible when is finite as long as finite size effects 
are not too large. One would expect that the tunneling which restores symmetry at finite A^ 
to be of order exp(— A^^ • const.) which could be very small even for rather small values of A^. 
It would be very interesting to check this idea using lattice Monte Carlo simulations of these 
systems. 

Extremal corners of the phase diagram 

Finally we discuss extremal corners of the phase diagram, which are labeled 1,..,4 in Figure ^. 
Table 1 gives the values of k/7, A/7 and the values of the contributions to the energy density. 
Obviously the four extremal cases can easily be understood by the magnitude of the energy 
densities Pfiux, PAd and pp^p- Point 1 is in the extremal corner of the low density phase. All 
three densities are rather small. Points 2 and 3 are both in the high density phase, in areas 
which are dominated by adjoint charges (Point 2) and fundamental charges (Point 3). It is 
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nice to see, how the energy density of the sources is dominated by the contributions of the 
adjoint charges and fundamental charges respectively. Finally Point 4 is in a region where 
both the density of the adjoint charges as well as the fundamental charges is high and of the 
same magnitude. 





1 


2 


3 


4 




0.00025 


0.00025 


0.25 


0.25 


A/7 


0.000252 


0.25 


0.000251 


0.25 


PFIux 


1.033 X 10-6 


0.063260 


0.072505 


0.105058 


PAd 


4.064 X 10-9 


0.160997 


0.000170 


0.189518 


Pf,f 


2.008 X 10-6 


0.000401 


0.411543 


0.435337 



Table 1: Values of n/'j, X/'-f and the contributions to the energy 
densities at four characteristic points in the phase diagram. 



It is interesting to notice, that when pushing the Points 2,3 and 4 to even higher values 
of k/j and A/7, one finds that the contribution of the charges to the energy density grows 
faster than the contribution of the flux. These extremal areas correspond to /i ~ — 1 + e. 
Expanding the necessary condition (|3.18| ) and the densities (|4.6| ) - ( ^4.8| ) for ~ — 1 + e (use 
( [OT| )) we obtain 



pAd 



A 

7 



1 + 




1 + o(yV7 + V7) 



in the limit K/7, A/7 ^ 1. This shows that at very high energy densities, the contributions 
from the charges are dominating. 



4.4 Summary 

In this paper we analyzed the thermodynamic properties of a model of static sources on a line 
interacting through non-Abelian forces. It was shown that the partition function takes the 
form of the partition function of a gauged principal chiral model. Using the eigenvalue density 
as collective field variable the Hamiltonian for the eigenvalue density in the large A^-limit was 
computed. We gave a static solution of the corresponding Hamilton equations. For the special 
case of only two types of charges, the static solution was parametrized using the parameter /z 
proportional to the Fermi energy. In particular the case of two types of charges transforming 
under the adjoint, and charges transforming under the fundamental representation of the 
gauge group was considered. Expanding the parametrized solutions at /i = 1 we established 
the existence of a straight line in the phase diagram where the free energy exhibits a third 
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order phase transition. We proved that the third order behaviour terminates at a second order 
point. The critical Une then continues as a first order fine, which was determined numerically. 

The whole phase diagram was interpreted by analyzing the contributions of charges and 
flux to the energy density. We found that for > 1 the system is characterized by low energy 
densities, while for /j, < 1 the densities are high. The behaviour of higher winding Polyakov 
loop operators provides a powerful tool for a further characterization of the high and low 
density phases. In the low density phase we found exponential suppression with increasing 
winding number, while for the low density phase we proved a power-law behaviour. 

It is tempting to analyze if the change of the behaviour of the Polyakov loops can be 
seen also in deconfining phase transitions in higher dimensions, and for finite N. If the 
picture persists, the Polyakov loop operators would provide a powerful criterion to analyze 
deconfinement phase transitions. In particular in lattice calculations this operator would be 
easier to analyze than other criteria as e.g. the area law of the Wilson loop operator. 

Also the 2-dimensional model could be generalized in several directions. It would be in- 
teresting to analyze non-static solutions of the Hamilton equations and different boundary 
conditions which might be used to include a ^^-term. Loop expansion of the fermion determi- 
nant of QCD2 with large quark masses could be used to relate the fugacities of the non-Abelian 
gas analyzed in this article to the mass parameters of QCD2. 
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Appendix 



A.l Cubic term in the Jevicki-Sakita Hamiltonian 

In this section of the appendix we bring the second term on the left-hand-side of Equation 
(|2.3CI| ) to a convenient form. The variation by p{9) of this term is given by (up to a factor) 



y\;{e) = [v j de'p{e') cot ^-y^ j -2v j de'p{e')de"p{e") cot ^-y^ cot ^' ^" 

In the integrals, we change variables to the complex variable 

t = e'S , t' = e''' , t" = e''", 
so that the integrals are over an interval on the unit circle and 

V J d9'pi9') cot = vj YT^P^''^ = ^ + 2^/ "^''jhp^''^ ■ 

We obtain 



w{t) =a(v [dt'^X-SV (dt'dfP^J^ 
^' \ J t-t'J J t-t't'-t" 



t'{t' - 1") 



The last two terms are constants (which because of the normalization condition the the 
density p{t) must satisfy, must be irrelevant), and the first term can be found by the following 
argument. We consider the function 



G{z) = Jdt 



Pit) 
t-z 



This function is analytic everywhere in the complex plane except on the arc of the unit circle 
where the eigenvalue density has support. Obviously 

G{z) — s> as \z\ oo . 

Also, by letting z approach the support of pit) from outside and inside the unit circle, we 
obtain (e > 0) 

liniG(t(l+e)) ^ G+(t) =vf dt'^ - mp{t) , (A.l) 

e^O J t — t 

and 

\imG{t{l-e)) = G4t) = V J dt'p^^ + nxp{t) , (A.2) 
respectively. The function 

pit) ^ f Pit') 



K{z) = G\z) - 2 J dt-^V J dt'^ 



t 
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is obviously analytic everywhere except eventually on the support of p. Using ( [A. 1| ) and ( [A.2| ) 
one finds that it is continuous across the support of p since 

ImiKitil+e)) - \imKitil-e)) = 0. 

Thus K is an entire function of z. Furthermore, since it vanishes at infinity, and is analytic 
everywhere 

K{z) = 0. 

Then, taking the real part of K on the support of p gives 

W(t) = 47rV^(t) + const. . 

This is the functional derivative of the term 

47r2 



J dOp^iO) + const. , (A.3) 
which is proportional to the second term in the Hamiltonian in the Schrodinger equation 

<m). 



A. 2 Analysis of the parametric integrals I{fi),M{fi) 

Here we analyze the properties of the parameter integrals I{p) and M[p) originally defined 
in (|3.15|) . A transformation of the integration variable brings the integrals to the form 

4 



lip) = -f Jp + cos(9) B.(p + cos(e)) de , 

77 JO ^ 

M{p) = - r coseJp + cos{9) Rip + cos{9)) d9 , (A.4) 

TT JO * 



which shows, that they are independent of J. H(..) denotes the step function. For p = I they 
can be evaluated explicitly, giving 

/(I) = ?^ , M(l) = 15^ . (A.5) 

vr 6 IT 

For p > 1, the range of integration is [0,7r], while for p < 1, the step function reduces the 
range of integration to [0, arccos(— //)]. I{p) can be expressed in terms of complete E{p) and 
incomplete E{p; 0) elliptic integrals 

^ J sySHf; (^31) tor ,< > 1 

= I »VmE (^31; aJin (;i±r)) for ^ < 1 ' 

When expanding I{p) and M{p) 

p = l±e,e>0 , I{l±e) = +6±I , M{1 ± e) = M{1) + S±M , (A.7) 
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one can use for /(/x) the well known formulas for the expansion of elliptic integrals (see e.g. 

6±I 



p3) to obtain 



Te\n{e)\ ± £-^[5 ln(2) ± 1] + o{e) . 



The variation 5±M can be related to 5±I by using 

-/ 



M(^) + 



[/i + cose]^/2 H{^i + cos{9))d9 . 



Expanding the left hand side using ( [A.7|) and the right hand side using Taylor expansion, one 
obtains 

— = —3—— ± -e . (A. 9) 

M(l) /(I) 2 ^ ^ 

We remark, that and M(/i) have first derivatives that diverge logarithmically as fi ap- 
proaches 1. This fact can be seen from (|A.6| ), (|A.9|) and known formulas for the derivatives 
of elliptic integrals. 

Finally we denote the large fi behaviour and the behaviour at ~ — 1. In both cases the 
results can be computed easily by expanding ( [A.4| ). We obtain 



M(/i) 



v^4[l-ll + 0(l//i^) 



(A.IO) 



for 1, and 



n-i+e) 

M{-l + e) 



eV2 + o{e) 
eV2 + o(e) 



{A.ll] 



for < e < 1. 



A. 3 Large k behaviour of Mfc(/i) 

In this section of the appendix we analyze the large k behaviour of the higher moments Mk{fi) 
given by 

Mkifi) = - r cos ke^fi + cos{e) H(/i + cos{e)) dO . (A.12) 
Again it is straightforward to evaluate the integrals at /i = 1. We obtain 

A4(l) = i^(-l)'=+'^ . (A.13) 

For /X 7^ 1 no exphcit solutions are available. For the case > 1 we analyze the complex 
contour integral 

Ic = - I z''-\L+hz + z-^)dz . (A.14) 
m Jc V 2 
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Figure 5: Contour for the evaluation of the integral (A.I4) 



For > 1, the argument of the square root 

can become real and negative, if and only if z is real and negative, i.e. z = —r, r > 0. In 
particular f{—r) is zero ior r = r± = fi ± ^J [i^ — 1, and f[—r) is positive for r_ < r < r_|_. 

On the complex plane cut along the negative real axis \J f{z) exists uniquely and is holo- 
morphic. If C is a closed curve not crossing the cut, then Ic vanishes due to Cauchy's theorem. 
We choose a curve as shown in Figure ^. The outer circle is the unit circle, the inner one has 
radius ^ 

r(/i) = + . 

The two pieces parallel to the real axis {z = x + is) cancel in the limit e ^ (Note that 
f{—r) has zeros at r = r±. r(/i) was chosen to lie half way between -1 and — r_, so that both 
integrals can be expanded in e showing that they cancel each other.). 

Thus we find that the integral along the inner circle equals the integral along the unit 
circle (both integrals evaluated in counterclockwise direction). Mfc(/i) is given by the real 
part of Iq evaluated along the unit circle. Using the fact that the integral along the inner 
circle is bounded, we establish 

\MM\ < (r(/i))'i=^2/i + r(/x)+r(/i)-i. (A.15) 

This proves that for > 1 the moments Mk{fi) are exponentially suppressed with increasing 
k. 
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Finally we notice that for — 1 < /i < 1 exponential suppression is excluded by well known 
results of Fourier transformation. Using the fact that Mk{fi) is proportional to the Fourier 
coefficients of \f + cos(^) we find 

1,. / X , 



/x + cos(^) = -Mo{fi) + Y.^^cos{ke) 
^ k=i ^ 



The derivative with respect to Q then is given by 



^ ^/i + cos(^) 



fc=l 



If the Mk{fJ^) were exponentially suppressed, the series on the right hand side would converge 
to a finite constant for all 6. This contradicts the fact, that the left hand side diverges at 
9 = arccos(— /i) for /i < 1. Thus for — 1 < < 1 (for /i = 1 see the result (|A.13| )) we find only 
power-like suppression of Mfc(/i) for large k. 
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